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ABSTRACT 

We investigate a one-dimensional loop model for the evaporative 
cooling of the coronal flare plasma. The important assumptions are 
that conductive losses dominate radiative cooling and that the evapora- 
tive velocities are small compared to the sound speed. We calculate 
the profile and evolution of the temperature and verify that our assump- 
tions are accurate for plasma parameters typical of flare regions. 

The model is in agreement with soft x-ray observations on the evolution 
of flare temperatures and emission measures. The effect of evaporation 
Lc to greatly reduce the conductive heat flux into the chromosphere 
and to enhance the EUV emission from the coronal flare plasma. 
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I . Introduction 


It is becoming widely accepted that flares evaporate large masses 
of chromospheric material into the corona. The most compelling arguments 
for the occurrence of this process come from soft x-ray observations. 
Measurements of the emission measure, which appear to be well over 

10^ cm for some flares (Phillips and Neupert, 1973)» indicate a 

T 16 

total mass of hot (T ^ 10 K) plasma in the corona as large as 10 g. 

Similar mass estimates are obtained from loop prominence systems which 

are known to occur in post-flare regions (Jefferies and Orrall, 1965 ^. 

This amount of plasma is too large to condense from the corona and must 

somehow be extracted from the chromosphere. 

Neupert (I968) was perhaps the first to realize that corona flare 
plasma is the result of evaporation from the chromosphere. Similar 
ideas were expressed subsequently by Hudson and Ohki (1972) and by 
Sturrock (1973 )• In these papers, it was implicitly assumed that 
evaporation occurs during the initial deposition of energy into the 
chromosphere by means of high-energy electron streams. However, there 
appears to be observational evidence for evaporation during the decay 
phase of flares. The emission measure is often seen to increase while 
the total thermal energy does not increase and the temperature decreases 
(Kahler et al., I 9 TO; Horan, 1971). "Moving boundary" types of models 
have been prcooseJ, to explain these results, by Strauss and Papagiannis 
( 1971 ) and by Zauman and Acton (1974^. These models involve the assumption 
that flare heating occurs only at the top of a loop of coronal plasma, 
which results in an expanding hot region. The loop is assumed to be 
essentially infinite in length so that the effect of the heating on 

C. 



the chromosphere is not included. However, these models encounter 
certain difficulties: in particular, they cannot account for the origin 
of the large masses required by the soft x-ray observations. Further- 
more, if the "moving boundary" models were correct, flares should heat 
the corona well before they heat the chromosphere, whereas observations 
indicate that the Ho emission of flares begins well before the scft 
x-ray maximum. We therefore take the view that the increase of emission 
measure during the decay phase of a flare indicates that evaporation 
occurs not only during the impulsive phase of a flare but also during 
the decay phase. 

From the above considerations, it appears that there may be two 
evaporation mechanisms occurring during a solar flare. One if. associated 
with the direct heating of the chromosphere by the energy released 
during the impulsive phase of a flare, the energy probably being trans- 
ferred by an electron stream. This stage is likely to involve large 
pressure gradients and large velocities, quite possibly leading to 
shock waves. It seems likely that a separate evaporation process may 
occur after the primary energy release is finished, this evaporation 
being driven by the large conductive heat flux from the high-tempera- 
ture flare plasma contained in magnetic flux tubes above the chromosphere. 
In this stage, evaporation is a cooling process competing with radiation. 
It seems likely that energy loss by evaporation dominates early in the 
decay phase, and energy loss by radiation dominates later, when the 
density of flare plasma is higher and the temperature is lower. More- 
over, it seems likely that this evaporation stage would involve only 
small pressure gradients and velocities. It is this later, "gentle", 
evaporation process which we discuss in this article. 
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II. Model 


Since the magnetic field strength in typical post-flare regions 
appears to be large (B ^ 10 g. Rust and Bar, 1973 ), we assume that it 
dominates the plasma so that mass flow and heat transfer are both 
parallel to the magnetic field. The geometry which we adopt is that 
of a flare loop consisting of a narrow current-free flux tube with the 
independent variable s measuring distance along the loop (parallel to B) 
from the top of the loop. The cross-sectional area of the flux tube 
is given by A(s) and will depend on the form of the magnetic field. 

The field is assumed to be that of a point-dipole or of a line-dipole, 
but we also calculate the case of a loop of constant cross section 
corresponding approximately to the familiar plane-parallel model. 

We make two assumptions on the state of the plasma. The first is 
that radiation losses are negligible. We later check our assumption 
by comparing the radiative cooling time with the conductive cooling 
time to be derived below. 

Our second assumption is that the velocities generated by the 
evaporation process are small compared with the speed of sound. Since 
the temperature is high enough that gravitational effects are negligible 
(i.e., the scale height is large compared with the size of the flux 
tube), this assumption is equivalent to setting 

If-O . (2.1) 

. -i.iL. it Sfcfcni paradoxical that we can calculate the 

plasma velocities by neglecting the pressure gradients, since these 
gradients create the velocities in the first place. However, one 


h 



finds that, if the gravitational term is negligible, the fractional 
variation of pressure along the flux tube is of order M^, where the 
Mach number M is given by 

M = v/C , (2.2) 


where C is the speed of sound. Hence our assumption (2.l'i is equivalent 
to the assumption that M « 1, which must be verified a posteiiori . 

Equation (2.1) shows that the pressure is independent of s . We 
assume that there is no heat input into the loop, and we ignore 
radiation. Hence the total energy in the loop is constant, which 
implies that the pressure is also independent of time- 


p(s,t' = const. 


(2.3) 


On using equation (2.2), one finds that the energy equation may be 
expressed as 


. 1 ^ = i ^ (aK 

2 p \Bt 9s / A 9s \ 9s / 


(2.4) 


On combining equation (2.4) with the continuity equation 


I? + i I; * 0 • 


(2.5) 


we obtain 


h p(|p'' ■ ■'If)] - ° 

We see from (2.6) that 

a(|pv - = f(t) . (2.7) 
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The temperature of the chromosphere (the base of our model) is 
sufficiently low that the heat flux F, given by 




F = - 



( 2 . 8 ) 


may be neglected. Similarly, the density there is so high that the 
mass velocity v may be neglected. These boundary conditions thus lead 
to the simple relation 


_ 2 K ^ 

^ 5 P Bs ’ 

which may alternatively be expressed as 



(2.9) 


( 2 . 10 ) 


Equation (2.10) can also be obtained from the symmetry of the loop since 
if the loop is symmetric about the top, it follows that F and v must 
vanish there. 

The heat equation may now bo expressed with only temperature as 
dependent variable. On replacing v in equation (2.4) by the expression 
( 2 . 9 ), and then eliminating p in favor of T by means of the equation 
of state, we obtain 


We use the Spit’.er 


5. £ ^ ^ i L- 

2 T Bt A 9s 



( 1962 ) formula for thermal conductivity, 


(2. IP 


K 


a T 


5/2 


( 2 . 12 ) 


where, with sufficient accuracy for our purposes, we may take (V = 10 
Equation (2.11) is nonlinear, but it is amenable to solution by 
separation of variables. We assume that T is of the form 
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(2.13) 


T(s,t) = T^e(t)cp(s) . 

For convenience, we choose 


0(0) = 1 , cp(0) = 1 


(2.1U) 


so that T, is the maximum temperature of the plasma (the initial 
M 

temperature at the top of the loop). We may substitute (2.13) into 
equation (2.11) and obtain 


2 


otX 


72 


^ ' A dl (''d^) ■ 7 ♦ 


M 


where >ll(s) is related to cp(s) by 


^ = 



(2.16) 


In equation (2.15). i^he variables are separated so that we may 

2 

set each side equal to -k. (independent of both t and s). The 
equation for 6(t), 




(2.17) 


may be integrated to give 

0(t) = (1 + , 


(2.1R) 


where 


T = 

2oT 

The equation for If becomes 

d^'ll 2 1 . 1 dA dji 

j_2 ■ 7 i;i Us) A ds ds 



(2.19) 


( 2 . 20 ) 


7 



We see from equation (2.14) and the symmetry of the flare loop that 
we require that 

\(l = 1, ^ = 0 at s = 0 . (2.21) 

However, if is the value of s at the "base” of the model (i.e., at 
the chromosphere), the temperature must drop to the low chromospheric 
value at s^^. For our purposes, it is convenient to set 

m(sj^) = 0 . (2.22) 

The position of the chromosphere will actually be not at the point 

4 

where il( = 0 (i.e. T = O), but where T s» 10 K, but this difference 
is small compared with the size of typical flare loops. The conditions 
(2.21) and (2.22) can be met only for a unique value of k in equation 
(2.20\ This value will of course depend implicitly also on the area 
function Afs\ For either type of magnetic-field tube (that produced 
by a point dipole or that produced by a line dipole), we shall wish 
to consider a range of values of s^^. Rather than prescribe S|^ and 
then determine k as an eigenvalue of the differential equation (2.20 ), 
it is clearly simpler to consider a range of values of k and determine 
Sj^ for each value. This may be done by finding a solution of the 
equation with the initial conditions (2.21), and then finding the point 
s at which \1( = 0. 
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III. Results 


(a) Temperature Profile 

For the case of a loop of constant cross-section, equation (2.20 ) 
can be integrated directly yielding: 


XTTTi . 750) 

y3 + 1 - vAo 3 


( 3 . 1 ) 


+ 2.3^^^ e(3(w), sin 75°) , 


3 (cp) = arc cos 


( 3 .?) 


where 

x /^ - 1 + \/tt 

,\/3 + 1 - y<p 

and F and E are the elliptic integrals of the first and second kind 
respectively. The constant k can now be expressed in terms of the 
position Sj^ of the base. Using equations (3*1) and (3.2), we find that 
for the planar case, 


k = 1 . 60 /s, 


(3.3) 


For loops of variable cross-section, equation (2.20) has been 
integrated numerically. The important parameters that characterize the 
geometry are the loop length Sj^ and the ratio of the cross-sectional 
area A(o) at the top of the loop to the area •^(sj^) at the base. We 
term this ratio the "compression factor" F: 


r = A(0)/A(s^) = B(Sj^)/B(0) . ^3.^-) 


Rust and Bar (1973) find that F ^ 30 in the flare of August , 19’^2. 


ccfs ) . 


For a given form of A(\s), we may solve equation ^2.20) for 
In Figure 3.1 we plot the temperature scale height H implied by our 
model , 


H I'T ) 


1 

T ds I 


( 3 . 5 ) 


as a function of temperature along the loop. This is done for the 
case of a constant-area loop and for a flux tube due to a line-dipole 
source with a compression factor T = 30. The results for a point-dipole 
source are almost identical to those of a line-dipole. The fact that 
the scale height H becomes infinite at the top of the loop simply 
reflects the fact that the spatial derivative of the temperature is 
zero at the top of the loop. This peculiarity is without physical 
significance; it simply means that the spatial var..ation of temperature 
is determined by the second derivative of the temperature rather than 
the. first derivative. 

The true size scale over which T varies is, of course, less than 

the loop length. For T < 10 \ the scale height is small compared to 

”2 5/2 

the loop length, H(T) < 10 s^^, and decreases as T"^ in both the 

constant-area and the variable-area cases. The loop would therefore 
appear to an observer to be almost isothermal with a very thin transi- 
tion region, as is the case for the quiet corona. 

For comparison we also plot in Figure 3.1 the temperature scale 
height given by the static model discussed by Antiochos an»- oturrock 

( 1 07(3 ) , From this figure we find that at lew temperatures (T < 1 

2^5 

the temperature profile is given by: T s for evaporative cooling, 
and T « s for static cooling. Hence the heat flux F is linearly 



proportional to T in the fvaoorative models, whereas F ~ const, in the 
static models. This result argues strongly in favor of the evaporative 
models . 

For temperatures and size scales typical of flares, the heat flux 
into the chromosphere implied by the static models is very large, in 

3 -2-1 

excess of 10 ergs cm s . It is difficult for static models to 
account for the dissipation of such a large flux; however, in our 
present model this flux is reduced by the evaporation process. The 
heat flux into the chromosphere is approximately three orders of 
magnitude lower for an evaporative model than for a static model. 


(b) Cooling Time Scale 

We can evaluate the cooling time in equation (2.13) by using 
equation (3.3)- For typical loop parameters, 

I ~ 10'^ K, n ~ 10^^'^ cm"^ and s^^ ~ 10"*^ cm , (3.6) 

equations ( 2 . 19 ) (3*3) yield 


2 

T ft; 10 sec . 
c 




The radiative cooling time is given by 


T 


R 


1 2 

^ n^A(T) 


(I.')) 


where A(T) has been computed by Raymond et al. (19"'6). For the above 

o 

parameters (3.6), we find that T 10"' sec; hence, our assumption that 

K 

conauction dominates the cooling is justified. However, as the tempera- 
ture decreases the ratio t /t„ increases until eventually radiation 

dominates and evaporation becomes unimportant. For the loop described 

r 

above, we expect this to occur at T < 5 ^ K. 
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We .lave also calculated the cooling times for loops due to either 

a point dipole or line dipole source. The effect of the geometry is 

shown in Figure where we plot the ratio of the cooling times for 

a 1 ''op and for a planar (constant cross-section) model as a function 

of r. As in the static case (Antiochos and Sturrock, 1976), the effect 

of the compression factor is to increase the cooling time; however, 

we still expect that conduction dominates over radiation during the 

7 

initial stages if T ^ 10 K. 

The other main assumption in our model is that the evaporation 
velocities are small compared to the speed of sound. For the planar 
case, the Mach number can be obtained directly from equations (2.9) 
and (3*1 ■ : 

c 
2/3 

The maximum value of M occurs at t = 0 and cfl = (2/5) .5« For 

the parameters in (3-6), we find that M ^ 10 \ so that our assumption 
is valid. From the form of the time-dependence of M and T, we find 
that M T^ for fixed s, hence our aporoximation improves as the 
loop cools. 


(1 + (3.9) 


= - = lo' 

c 


■ 10.6 




1/2 


C5 (1 


ps 


( c ) Rad i a t i on 

Various authors have attempted to calculate the emission measure 
and the temnerature of flare regions from soft x-ray observations 
^e.g. Kahler et al., 19T0; Horan, 1971^. Craig and Brown (1976) 
point out that these calculations are inherently inaccurate, and 
suggest that one should instead assume a definite model for the 
structure of a flare region and then test this model by comparing the 
expected radiation with observations. We, therefore, calculate the 




tmA v; A T , 


1 ? 


evolution of the soft x-ray emission predicted by our evaporative loop 
mode . 


Generally, the observations consist of measurement of the total 
flux of x-rays in some energy range. Hence we must compute 


l(E^, E^, t) 



ri(E, X, 


t I 


( 3 . 10 ) 


where the spatial integration is over the volume of the loop and r) 

is the energy flux at the earth from a unit volume of plasma. If the 

dominant radiation process is bremsstrahlung, the emissivity ri(keV 
-3 -1 -1 

cm s keV ) is given by 


T| = 


^ -E/kT 2_,-l/2 
ae n 1 


(3.11 ) 


where a = 
(3.1) and 
a loop of 


10 for T in °K, (Kahler et al., 1970). Using equations 
( 3 . 11 ) and integrating equation ( 3 .IO) over E, we obtain for 
constant cross-section A: 


X(r^, r^, t) 



(3.12) 


where 


and 


= 10 


■ 13.8 

t,3/2 

Si 


^ = 


"1 


"2 


(3.13) 


(3.14) 


In Figure 3.3 we plot L for the particular values rj^ = 1/2 and 
1 as a function of ^/t^. Assuming a maximum loop temperature 
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= 2 X 10 K, these values for and correspond to a soft x-ray 
range from 1-2 keV. We also plot the time dependence of the temperature 
and of the eiission measure in Figure j.3» Note that even though the 
temperature is continuously decreasing, Z increases initially due to the 


increasing emission measure and has a maximum at = 50 . For larger 

values of r^ and r^, e.g. r^ = 2 (E^ = 4 keV) and r^ = 4 (E^ = B keV ' , 
the exponential term in the emissivity (3.11) dominates, and L decreases 


monotonically . We find that these results are insensitive to che 


compression factor T because the x-ray emission comes mainly from the 


high-temperature material high in the loop. 

The qualitative features of Figure 3.3 are in agreement with soft 
x-ray observations. However, we would not expect to find detailed 
agreement, since we have calculated the emission from only a single 
loop, whereas a flare is likely to consist of many loops at different 
stages in their evolution. In addition, we have neglected any heating 
which may be present during the decay phase of flares. However, our 
results do predict a rise in soft x-ray flux concurrent with a decrease 
in temperature, whereas in the static models the emission must decrease 


with the temperature. 


We believe that a strong observational test of our model may lie 


in the spectrum of the radiation. From the form of the scale height. 
Figure 3. It we note that an evaporative loop contains relatively more 
plasma at low temperatures 10 K) than does a static one, which 
Implies that the spectra of the emitted radiation (especially at EUV 
wavelengths 1 should be significantly different in the two cases. I’his 
is clearly illustrated in Figure 3.^*^ where we plot, for both models. 
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the total radiative losses e (X ) due to plasma at different temperatures 
along the loop, i.e. 

€(I) = n^A(T) H(T) , (3.15) 

where we use Cox and Tucker's (I969) values for A(T^, and H and n are 
in units of loop length s^^, and density at the top of the loop n"s = O), 

respectively. The temperature at the top of the loop is assumed to be, 

7 

T,, = 10 K. We note that in the static case the highest temperature 
M 

plasma is responsible for most of the energy loss whereas, in the 

5 

evaporative loop, plasma at ~ 5 ^ 10 K radiates most strongly. In a 
later paper we intend to calculate the spectra for both models in detail, 
in particular the emission in various EUV and x-ray lines, and to compare 
these calculations with recent Skylab observations. 
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Figure Captions 


Figure 3. 


Figure 3. 


Figure 3. 


Figure 3. 


: The temperature scale height as a function of temperature 

along the loop for three models. The top solid line refers 
to an evaporative model of constant cross-section. The 
broken line refers to an evaporative model for a loop due 
to a line-dipole source with F = 30. The bottom solid 
line refers to a static model of constant cross-section. 

: Ratio of cooling times for a variable cross-section model 

as a function of F. The solid line refers to a loop due to 
a point-dipole source and the broken line to a loop due to 
a line-dipole source. 

: The ev »lution of the soft x-ray emission £ between 1-2 

keV, the temperature T, and the emission measure (E.M. ) 

7 

for a loop of initial temperature, T^^ = 2 X 10 K. 

: The radiative losses per unit logarithmic temperature interval 

of plasma in a loop as a function of temperature along the 
loop. The top line refers to an evaporative model and the 
bottom line refers to a static model. 
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Figure 3*1 
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Figure 3*3 
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